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3. SUMMARY 


An Ultrasonic Simulation of Orbiter, and it's Microwave Scanning Beam 
Landing System (MSBLS) , the runway and rain environment were performed in 
order to determine the multipath fading from the ground, and its possible 
degrading effect on the orbiter received beam shape and the associated land- 
ing guidance parameters. 

The on-shuttle antenna pattern of the MSBLS receiver, -40® ^ <{> ^+40°, 
-23° ^ 0 ^ +27° , and the azimuth and elevation beamwidths 7° EL x 2° AZ and 
1.3° EL X 40° AZ respectively, were simulated by their corresponding ultra- 
sonic transducer beams. The scanning rate for the azimuth and elevation 
beamswas 1.75 degrees /second. The results were adjusted for full-scale 
maximum sinusoidal scan rates of 691° and 377° /sec for AZ and EL respective- 
ly. The rain drops were simulated by air bubbles, with a similar size 
distribution, in water. The rain volume was created along a part of the 
propagation path, and not on the runway, because it was found difficult to 
avoid an accumulation of bubbles on the runway surface and surroundings 
simulated by the model surface. 

The following are a summary of the results of the work reported in this 
final report : 

1. No-rain multipath was found to occur in front of the azimuth and 
elevation transmitters, as was the case in NASA aircraft flight 
tests at Edwards AFB in California, 

2. The ground multipath seen by the actual 50° EL x 80° AZ shuttle 
receiver coverage antenna is expected to be larger in amplitude 
than that measured by its simulated narrow beam antenna used 
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aboard the test aircraft in flight tests at Edwards APB, California, 
based on ultrasonic simulated results. 


3. All the multipath seems to occur at points on or in the immediate 
vicinity of the runway, and therefore the surrounding terrain does 
not contribute to multipath conditions for ranges up to 20,000 ft. 

4. Since the multipath reflections from the runway and its Immediate 
vicinity occur at angles less than or around 16 degrees from the 
azimuth, the reflection coefficient for vertical and horizontal 
polarization are within 3 db for this range of grazing angles and 
therefore ultrasonic simulation- of MSBLS is valid. 

5. The ultrasonically simulated multipath causes a -4 db beam center 
shift (caused by the very presence of runway as opposed to free 
space propagation between Orbiter and Azimuth scan transmitter) 
ranging from a minimum of 0,083 to a maximum of 0.625 degrees when 
referred to the free-space condition case. A similar shift is ex- 
pected for the case of elevation scan. 

6. The presence of simulated severe rain causes a further -4 db beam 
center shift by up to a maximum of 0.08° at 27,250 ft. range for 
azimuth scan, and 0.13° at 32,000 ft. range for elevation scan. 

The angular positioning accuracy of WPL - UH system is ± 0.01,^ 
whereas range and altitude settings were within ± 8 feet of the 
full scale value, 

7 . Chase aircraft in different positions around the shuttle seemed to 
produce no degradation of the received beam shape. 
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It is further recommended that; 


a. a cross correlation between full-scale beam center shift data and 
data obtained by ultrasonic simulation be obtained, 

b. rain effects on the Orblter received beamwidth be measured and 
analyzed in a test at Kennedy Space Center with MSBLS, so that rain- 
caused errors, if any, are clearly separated from the no-rain multi- 
path conditions and 

c. that cross-correlation be calculated between rain-caused multipath 
effects obtained for full-scale and ultrasonically simulated cases. 
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4. MSBLS DESCRIPTION AND IT* S ULTRASONIC /MODEL 


The Orbiter MSBLS system modeling consisted of three separate tasks 
namely the azimuth and elevation scan transmitter as well as Orbiter receiver 
beamwidths, the geometry, and the runway and associated terrain. Each of 
these areas are discussed below in addition to polarization and signal wave- 
form. 

BEAMWIDTHS OE TRANSMITTERS AND RECEIVER 

A summary of beamwidths of the azimuth and elevation scan transmitters 
as well as the Orbiter receiver are given in Table 4.1, and the corresponding 
beamwidth transducers are used in the ultrasonic model. One half inch dia- 
meter cylindrical disc transducers with specially constructed lenses were 
used to simulate these beamwidths. The Orbiter receiver coverage was simu- 
lated by 80° X 50° (AZxEL) transducer whose two diagonal pattems/cuts are 
shown in Fig. 4.1. This unit satisfies the 4 db/8° period ripple and beam- 
width requirements. This was the most difficult beam pattern to design, 
since procurement of an exact 80° x 50° pattern requirement would have been 
prohibitive in coat and time, and therefore this was considered to be the 
optimum from cost and delivery time standpoints. Since the Orbiter receiver 
coverage shows a considerable ripple therefore its transducer coverage re- 
presenting it is probably as good a fit as one can achieve in practice. 

The 40° X 1,3° (AZxEL) elevation scan transmitter was simulated by its 
corresponding transducer (XDCR) with both pattern cuts shown in Fig. 4.2, 
whereas the 2° x 7° (AZxEL) azimuth scan transmitter cpunterpart transducer 
patterns are shown in Fig. 4.3, A close examination of these two coverages 
shows that these approximated the full-scale antenna patterns. 
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Thus having approximated the beamwidth of the antennas, one comes to the 
scan rates. The maximum full-scale sinusoidal EL and AZ scan rates were ap- 
proximately 377® and 691° /second respectively, and such rates are well nigh 
impossible to duplicate in a water tank for a variety of reasons. Thus, 
practicality of slow scan rate, necessitated by the requirement of disturbing 
the propagation medium (water), of the order of 1.75° per second were used. 
The effect of slow scan rates is to preserve the high-frequency time varia- 
tions of the signal, whereas the high scan rates would tend to act as lowpass 
filter which filters out such high frequency content and smooths out the 
waveform distortion. This, when translated in terms of the Orbiter received 
beam-shape, means that any distortion of the said beam-shape will be smoothed 
out but the basic "dent" shall remain, or in other words frequencies of the 
order up to 215 and 395 times for EL and AZ respectively would be filtered. 
This point is also discussed in the chapter on rain simulation and discussion 
of results and therefore no further details are provided in this chapter. 

LANDING TRAJECTORY AND GEOMETRY OF TRANSMITTERS /RECEIVER 

During the first phase of this work, it was mutually decided to select 
points on the Orbiter landing trajectory provided by NASA, and to statically 
simulate the received signal beam-shape for analysis of the effects of the 
ground multipath. The location of the elevation and azimuth scanning trans- 
mitters and a summary of their beamwidths and orientations are shown in Fig. 
4.4. In each of these figures, model distances calculated using linear 
wavelength scaling are also shown (See Table 4.1 and Chapter 5). Linear 
wavelength scaling preserves the phase and amplitude effects exactly except 
for a fixed db shift between the electromagnetic and ultrasonic received 
signal levels. 
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In the case of the Orbiter trajectory, it was essential that both the 
elevation and range as well tilt angle of the receiver with respect to trans- 
mitters be maintained exactly the same as in the full scale electromagnetic 
case. Therefore, in order to simulate exactly the ground multipath, one 
uses a simple geometry shown in Fig. 4.4 in order to arrive at a range scale 
factor of 3200:1 and an elevation scale factor of 3310:1. Based on these 
scale factors, and the actual grazing angles of multipath, one calculates 
the corresponding range and elevation for the trajectory in the water for 
the ultrasonic model. The results of such a calculation are shown in Fig. 

4.5 where range, elevation, beam depression angle are listed for both the 
full-scale electromagnetic as well as for the ultrasonic model case. 

POLARIZATION MD SIGNAL WAVE FORM 

Since MSBLS uses vertical polarization and that the grazing angles of 
the order of or less than 16® are involved in the ground multipath mode, it 
is not unreasonable to expect a difference of two to three db in the reflec- 
tion coefficients for vertical and horizontal polarization when one combines 
the phase and amplitude effects of the surface roughness. As the ground 
roughness decreases or it becomes a better conductor, the differences be- 
tween amplitudes for two polarizations reduce to negligible proportions (see 
Chapter 5) for such small grazing angles. This is even much more valid for 
rain conditions, since water cannot support shear (lateral) waves, longi- 
tudinal (compressional) waves are used to simulate the vertical polarization. 
Further details on this matter are provided in Chapter 5, with examples. 

Continuous sine wave signals was transmitted and received in the ultra- 
sonic model system, and the other details of the experimental system are 
given in Chapter 6. 
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OTHER FACTORS AND RUNWAY MODEL 


Chase planes were simulated by l/72th scale commercially available jet 
aircraft coated with appropriate aluminum-zinc loaded epoxy in order to 
guarantee appropriate Fresnel reflection coefficient, 

The runway was constructed of lucite— air— lucite sandwich which was com- 
pletely air-sealed. The Fresnel reflection coefficient of the smooth surface 
was found to match that of paved runway or a lake bed so well in the 0-30° 
grazing angle range that this construction was preferred over our classical 
built-up models used successfully in the past decades. Photographs of the 
deck are shown in Chapter 6. 

The area' surrounding the runway was buffed to make it appear rough for 
2.25 megahertz ultrasonic waves used in the model experiment based on X/lO 
roughness criterion for vertical Incidence (Hayre 1962) . A general reflec- 
tion coefficient check of the model runway and the surrounding area were 
verified before the final tests were run in the model experiment. 

The following symbols are used in Table 4 following this page: 


L = Length 

f = frequency in air Cuctual frequency) 
f^ = frequency in water (model frequency) 

R = Range scale factor 

c^,c^ = velocity of wave propagation in air, in water 

1 X = wavelengths in air and water 
a w 


Subscript a = for air 

w = for water 
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TABLE 4.1 EDLLSCALE AND MODEL MSBLS PARAMETERS 


1(a) EULLSCALE ORBITER 

Length = L = 1467.06" 

Wing Span = 936.68 

Height Stab. Tip = 679" 

f = 15.4 - 15.7 GHZ 
a 

Center f =15.55 GHZ 
a 

Q 

C = 3 X 10 m/sec 
a 


ULTRASONIC MODEL 
50.76" 

32.4" 

23.5" 

Ultrasonic wave velocity 

3 

in water = c = 1.5 x 10 m/sec 
w 

Ultrasonic frequency f = ? 


WAVELENGTH SCALE FACTOR (R) 

R = X /X = (Ca/Cw) . (fw/fa) = 28.938 29 

a w 

L 'v 48" -J- 50" ^ R = 29 f = 2.25 MHZ 

w 

Wing Span 24" 36" 


1(b) MSBLS ANTENNA/RADAR SIMULATION 

Elevation Beam 1.3“ EL x 40“ AZ (-3 dh) (0-30“ EL. SCAN.) 

side lobe -18 db 

Azimuthal Beam 7“ EL x 2“ AZ (-3 db) ±- 15“ scan in AZ) 

side lobe -20 db 

SHUTTLE MSBLS RECEIVER ANTENNA 

-40“ 4 ^ < +40“ , -23“ <04 +27“ 

Ripple 4 db ma/8“ period 

Gain grad.: a maximum of 0.5 db /degree 

Ultrasonic Beams Identical 
Transducer 1/2"; 2,25 MHZ, Special Lenz 

1(c) Controlled Azimuth and Vertical Scan Using Flexible Shaft/Scientific 

Atlanta Positioner 

Position Accuracy - 0.01“ 

1(d) Polarization - Vertical (Longitudinal Waves) 
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EL C+27® 


FIG. 4,'l ORBITER HSBLS RUNWAY CONFIGURATION 






5. VERTICAL POLARIZATION MSBLS MtILTIPATH A3 nTO ITS 
ULTRASONIC SITUATION 


In case of acoustic waves In water, there is no polarization, and the 
pressure wave P is longitudinal. The associated particle velocity U is al- 
ways along the direction of propagation. It will be shown later that the 

pressure wave P is an exactly analog of the E type polarization. Since 

n 

liquids cannot support transverse waves, there is no parameter exactly anal- 
ogous to polarization in electromagnetics, except for airframe mounted 
antenna, near-grazing and near-vertical incidence cases, and for the case of 
almost perfectly conducting surfaces. 


WAVE EQUATION 

The wave equation for a uniform plane electromagnetic wave in free space 
and acoustic wave in water are analogous (Morse, 1968): 


E. M. Waves (Air) 


Acoustic Waves (Water) * 


n2ir 

V E = ye — ^ 
3t 


2 — 

V'^H 


ye 


3^H 

3t 


2 — 
V U 


= pK 


3^U 

3t^ 


V^P 


= pK 


3^P 

3t^ 


where; E = Electric field intensity vector, (volts /m) 
H = Magnetic field intensity vector, (amp/m) 

P = Scalar pressure field, (N/m ) 

U = Particle velocity vector, (m/sec) 

3 

p = Density of medium, (kg/m ) 

2 

k = Compressibility of medium, (m /N) 
e = Permittivity of medium, (f/m) 

» y = Permeability of medium, (h/m) 


(5-1) 


(5-2) 


5.1 


i wt 

Assuming e variation for both systems, solutions to the above cases 
of uniform plane waves propagating in negative z-direction are; (Kraus, 1973) 


Electromagnetic- 


— 'sr 1 (wt+kz) 

X 0 


y 0 


Acoustic 

P . p J ■ (5-3) 

z 0 

B .B.J (5-4) 

2C (J 


Since these forms of the solutions are identical, and differ only by the 
physical constants in the respective systems, therefore one concludes that 
there is an exact analogy between these two systems. 


BOUNDARY CONDITIONS 

Boundary conditions in electromagnetic case require that tangential com- 
ponents of E and H are continuous across the boundary; i.e. 


n X (E^ 

0 
It 

1 

1 

(5-5) 

n X (H^ 

- H^) = 0 

(5-6) 


For the acoustic case, the pressure P and normal component of particle velo- 
city II are also continuous across the interface between the t;wo media: (See 
Fig. 5.5) 


Pi = 1*2 " ^ 

n- (U^ - Uj^) = 0 (5-8) 

Boundary condition at a perfect conductor require the tangential electric 
field to be zero', and the tangential magnetic field to be two txmes the inci- 
dent field. On the other hand, at a perfectly elastic wall, the dynamic 
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acoustic pressure is zero and the normal component of the total particle 
velocity is twice the normal component of the incident particle velocity. 

In each of the situations there is no wave propagation beyond the interface. 

It has been well established, that specially constructed acoustic models 
with air-backed surfaces behave exactly like air /aluminum or air/steel inter- 
faces or air/good conducting ground or pavement for electromagnetic waves 
in so far as the Fresnel Reflection Coefficient are concerned. (See Fig. 

5-1) . This is also valid for NASA shuttle runway simulated by luclte-air- 
lucite sandwich. 


REFLECTION COEFFICIENT 

If a vertically polarized electric field vector is incident on a smooth 
surface in x-z plane at an angle 6^ along the vector, and is partially 
specularly reflected into the same dielectric medium and partially transmitted 
into the second medium at an angle ^2 shown in Figure 5-2, then reflection 
coefficient is, 

I t 

(5-9) 




^ 2-^1 
z; + z; 


whereas its transmission coefficient, T^ is, 


2Z„ Cos 0 
• ^2 1/2 

where: Cos 9, 


’ 1/2 
Z, = (-i)-^'^ Cos 9, 


(5-10) 


(5-11) 


(5-12) 
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In the case of horizontal polarization, the reflection and transmission 
coefficient for a plane interface between two infinite dielectric slabs are. 


(See Fig. 5.3) 


II II 




22-^1 


Z 2 + Z 1 


2Z, 

I-TT ~ II II 

Z^H-Z^ 


" ^2 1/2 

where; Z„ = (~) Sec 0, 
2 e„ 


" 1/2 

= Sec 9. 

1 J 


(5-13) 

(5-14) 

(5-15) 

(5-16) 


An acoustic wave incident in water on a solid surface splits into two parts; 
longitudinal and shear waves in the solid medium (see Figure 5-4) . The 
equivalent impedance of medium 1 and 2 are given (Brekhovsklkh, 1960) as, 


Impedance; Longitudinal, ~ *^i ^1 ®i 

Transverse, Z = p„ b„ Sec v„ 
t 2 2 2 


Longitudinal: Reflection 

coefficient, 


Z 

total 

^total 



Transmission 

p 2(Z cosZVj) 

coefficient, T- = — • -= — — 

^ ^2 ^total'^'^Zl 


Transverse Transmission 

p^ 2(Z^cosV2) 
coefficient , T , = — •= rz — 


(5-17) 

(5-18) 

(5-19) 

(5-20) 

(5-21) 
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where : 


(5-22) 



Z 


total 


= Shear velocity 
= Lame's Constant 
= Z2COS^2^2+2^.sin^2v2 


(5-23) 

(5-24) 


THEORY FOR SIMULATION OF VERTICAL AND HORIZONTAL POLARIZATION 

The Fresnel reflection coefficients for vertically and horizontally 
polarized electromagnetic waves at a boundary between two dielectric media 
are given as (Stratton, 1941) 




- U^k^cos 62 

+ y^k^Cos 02 


(5-25) 


^2^1^03 ©2 - u^k2Cos 

N' y2^2^o® ^2 y^k2Cos 0^ 


(5-26) 


where 0^ and 02 are angles of incidence and transmission respectively as shown 
in Fig. 5-2. The presence of conductivity in one of the medium modifies the 
reflection coefficient. 

If medium (1) is a perfect dielectric and medium (2) is conducting, then 
the propagation constants are defined as 

k^ = = w/y^e^ (5-27) 

^2 " “2 ^^2 

where a is the phase constant and 0 is the attenuation factor. 

For a dielectric /conducting interface, the Snell's law is modified as 
given below; 
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(5-29) 


Sin 0^ = k 2 Sin 0^ 

k a 

Sin 02 = k ° ••• 2 ~~ 2 ®1 (5-30) 

2 «2* h 

The above equation can be rewritten in a convenient form, given below: 


Sin 02 = (a-jb) Sin 0^ (5-31) 

The Cos 02 becomes a complex quantity as defined below: 

Cos 02 = A-(a^-b^-i 2ab) Sin^0^ Pe^^ (5-32) 

Then an application of the Snell's law yields: 

k2Cos 02 = A^-k^ Sin^02 '^,'^2~^1 (5-33) 


Now substituting Eq. (5-33) in Eq- (5-25) and Eq. (5-26), one obtains, 
(Stratton, 1941), the following: 


^1 = 


^ ia^-a^Qos e^)^ + $2 


^a2+a^Cos 6^)^ + 3^^ 


1/2 




( 02 - 62)^08 0^-aj^a2]^+l2a232Cos 


[(a^-S^JCos 0^4-0^02] ^+[ 2 a 232 Cos 


1/2 




(5-34) 


(5-35) 


In case of a conductor. 


2 _ ^2 T 

TI2 , 2 2 ^ 

w £2 


and 


«2 = 32 " 




(5-36) 


2 


(5-37) 



therefore and become, (Stratton, 1941), 



(1- Cos 0,) +1 

yj_«2 ^ 

^ l °^2 2 

(1+ Cos 0.) +1 

>'l“2 1 


|r^|^ = [x^+(2Cosej^-X2)^]/[xV(2Gos6^+x)^] 


(5-38) 


(5-39) 


Since in case of all conductors « 1, and denoting the quantity by 

h “2 \“2 

X, one can rewrite (5-39) as below; 



1 - 2x Cos 0^ 


(5-40) 


Similarly the reflection coefficient for vertical polarization, can be simpli— 


fied to the following form: 


2 2 
_ 2Cos 6- - 2xCos 0, + X 

Kl 

2Cos 0^ + 2xCos 0^ + X 

(5-41) 

At normal incidence (6^^ = 0°), jR^j = 1 - 2x 

(5-42) 

2 

, 1 „ 1 2 2 - 2x + X 

and |R^| = 

^ 2 + 2x + X 

(5-43) 

(1 + x) + 1 

(5-44) 


Applying binomial theorem to Eq. (5-39), jR^I reduces to 

l\l^ = l-2x - (5-45) 

which is the same as Eq. (5-42) for perpendicular polarization. For grazing 
incidence (0 = 90°), one obtains 
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(5-46) 


For near-vertical and near grazing angles, magnitudes of and are equal, 

and the range of angles in these two cases depends on the specific conducting 
material (depends on x) . Brekhovskikh , (1960) and Stratton, (1941) present 
plots of and from 0° to 90° angles for different interfaces. 

A sample calculation of R^. and R^ for air/aluminum interface (typical of 
air frame mounted antennas) for a range of incidence angles from 0° to 90° is 
given next. The value of x for non-magnetic materials is given by 


X 


^l”l _ ^1 

^2*^2 “2 


(5-47) 


The phase constant of E.M. wave in air, (a^ = for 

m = 10^ X 2 tt( 5-»-10) rps 
U = 1.25 X 10“® h/m 
= 8.84 X 10 “^^ f/m 
. *. = 104.35 ->■ 108.70 

The phase constant of E.M. wave in aluminum, 

9 

0 ) = 2tt X (5-»-10) X 10 rps 
V 2 = 1.25 X 10 “^ h/m 
O 2 = 3.5 X 10^ mho/m 
.'. U2 = 146.92 -s- 293.84 

• *^1 -4 

.'. X = — = 0.71 X 10 X 10 ^ 

“2 

Similarly for steel x = 0.122 ->• 0.174 
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efficient for the longitudinal waves from the model surfaces is also essen- 
tially constant near the near-normal and near-the-grazing angle domains for 
the angles of incidence encountered in MLS and in MSBLS multipath geometries. 

The basic angles of Incidence for MSBLS case for the range of altitudes 
considered here is approximately 70° to 82° , when measured from the outward 
normal. Since the sandy dry lake runway at Edwards may be said to have re- 
lative dielectric constant = 10 f/m, and the conductivity of approximately 
-3 

2 X 10 mhos per meter, thefefefe the lake bed and the metal runways both at 
Edwards AFB and Kennedy Space Center, under rain condition would somewhat 
behave like a material midway between a fresh water surface and the moist 
sandy surface. Since fresh water has an = 75, and its conductivity may 
vary from 1 to 20 mhos/m, it is reasonable to assume an average value of 10 
for the conductivity. 

Now, for an ideal smooth surface at 15 GHz, the Fresnel reflection co- 
efficients has been measured for a smooth sandy surface to be, for 70“ < 0 < 
82“ 


0,8 < < 0.91 

0.2 < < 0.4 

and for smooth fresh water surface to be 

0.83 < Sjj < 0.97 
0.5 ^ ^ ^ 0*1 

Under actual operating conditions, MSBLS microwave radiation would 

encounter a fairly rough surface in case of the dry lake bed, and, also under 
rainfall conditions in case of both the metal runways and the lake bed. 
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From extensive experience in scattering from random rough surfaces, it may 
be stated that 0.2 - 0.4 and 0.5 - 0.1 variations would considerably be 
smoothed out and one may expect the range of R^. to be within one to two dbs 
at the maximum. This condition is confirmed by the reflection coefficient 
variation for the model surface (shown in Fig. 5.1) of approximately 0.5 - 
0.75 db for 70® ^ 82 °, Therefore it is concluded that the special sur- 

face of water /air backed plexiglas (sandwich) interface simulates all MSBLS 
runway conditions reasonably well. 

For instance, most of the multipah from the runway comes from near- 
grazing incidence, and primarily the worst case is caused by strong specular 
reflections in the direction of the approaching shuttle, for which case 


l\l “ Ii>hI = l\l 


for rough surfaces 


(5-50) 


where >= longitudinal acoustic wave reflection coefficient. In fact even 
if js^l is not equal to |Ry.j and > as shown in Figs. 5.1 and 5.7, but its 
variation over the range of angles of Interest is relatively “insignificant, 
then the fading pattern at the ' receiver aboard the landing • craft is 
identical in shape to the full-scale electromagnetic case for either polar- 
ization. Of course the dynamic range of the fading shall be truly simulated 
by its ultrasonic model, and the amplitude of the simulated results would 
differ from the full scale by a constant factor, thus not affecting the -4 db 
beam centers used in this work. 

In summary it has been shown that |Fj^| simulates horizontal polarization 
exactly and the vertical polarization for near grazing and near normal angles 
of incidence to an'"accuracy of ± 1-2 percent. 
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FIG. 5,1 ACOUSTIC REFLECTION COEFFICIENT VS ANGLE OF INCIDENCE 
FOR H2O/AIR-BACKED PLEXIGLASS INTERFACE 
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FIG. 5.2 REFLECTION AND TRANSfllSSION OF VERTICALLY 
POLARIZED WAVE AT OBLIQUE INCIDENCE 



FIG. 5.3 REFLECTION AND TRANSMISSION OF HORIZONTALLY 
POLARIZED WAVE AT OBLIQUE INCIDENCE 
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FIG. 5.^ REFLECTING AND TRANSMISSION OF LONGITUDINAL 
ACOUSTIC WAVE AT SOLID/WATER INTERFACE 
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FIG, 5,5 BOUNDARY CONDITIONS AT AN INTERFACE FOR 
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FIG, 5,7 FRESNEL REFLECTION COEFFICIENT VS ANGLE OF INCIDENCE 
FOR AIR/WATER INTERFACE FOR (EHW). Rn aRy AND ACOUSTICS 


Horizontal (EMW) 


Acoustic Water/Air Interface 


INC I DENT ANGLE . DEGREES 






6. SIMULATION EXPERIMENTAL SETUP 


Since in the ultrasonic simulation, the propagation media is water 
volume, and therefore the simulation experimental setup is a mirror image of 
the full scale setup, in that the 4’ x 4’ model runway shown in Fig. 6.1a 
is turned over such that the highest elevation point on the trajectory is 
the deepest point in water. Two detailed views of the whole setup are shown 
in Fig. 6.1b and c. 

The ultrasonic frequency of 2.25 megahertz, as determined by linear 
wavelength scale modeling discussed in Chapter 4, is used to simulate the 15 
gigahertz signal used in MSBLS. Since the receiver coverage on Orbiter was 
used to design the receiver transducer, only the receiver transducer, instead 
of the Orbiter model with the transducer mounted on it, is deployed. 

A typical block diagram in Fig. 6.2 shows the electronic system. Either 
a 600 volt peak to peak Ahrenberg signal generator followed by an Ahrenberg 
tuner or a General Radio oscillator and HP 46 7A power amplifier are used to 
feed the azimuth and elevation scanning transmitters one at a time. The 
93 ohm coaxial cable connecting the oscillator to the transmit transducer is 
approximately 65 feet long, therefore it was found necessary to use the high 
voltage Ahrenberg source, in particular for the farthest (120") trajectory 
point. The receiver transducer is mounted on a special L-frame (See Figs. 

6.3 and 6.4) supported by cross and side braces in- order to minimize its 
displacement in water and to dampen its motion. The vertical leg of the 
L-frame is mounted on the side of a carriage, which can travel on two parall- 
el rails set up across the top of the 12 feet wide water tank. This two- 
rail system is mounted on an Otis elevator carriage which rides two rails 
mounted on the top of the two 42’ long side tank walls. This enables the 
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operator to adjust the horizontal range as well as the position, in azimuth, 
of the receiver unit. The mounting of the receiver unit is so supported 
that it can be raised or lowered to adjust the elevation of the receiver 
unit. Moreover, a special rotational axis system holding the transducer at 
the bottom end of this frame, is designed to be rotated from the top in order 
to adjust its tilt angle in the vertical plane. A protractor is mounted on 
the top of this mechanism for a direct indication of this tilt. This com- 
plex setup enables one to adjust all the necessary parameters of the receiver 
transducer. The operator has a list of settings of range, tilt angle, alti- 
tude (depth here) , and these are read on permanently mounted tape along the 
length of the tank protractor, and along the vertical leg of the special 
L-frame respectively. The accuracy of position is + 1/32", and the tilt 
angle for the transducer is ± 0.1 degrees. 

The 4' X 4* model runway is mounted, facing down into the water, on a 
specially designed frame, and coupled to a Scientific Atlanta positioning 
system mounted on top of 'the tank on a rigid platform supported by tank walls. 
The two scan units consist of plastic transducers in their respective hold- 
ers, mounted in specially constructed frames so as to provide the azimuth 
and elevation scanning capability. This is achieved by connecting a flexi- 
ble shaft from the positioning system shaft to the fixed shaft rotating 
transducer holder. Only one of the scan units was operated at any one time, 
whereas the other unit was removed in order to eliminate interference caused 
by the latter. The positioning accuracy of Scientific Atlanta system is 
± 0.01 degrees. 

The chase airplane models were hung by weight tied to a nylon fishing 
fiber below the model aircraft. Another tandem positioning connection was 
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provided outside the water in order to provide means to change their position 
relative to the shuttle receiver. These were manually positioned all over 
the space around the receiver, and were found to cause no multipath problems 
except when they were positioned in the common volume of the transmitter and 
the shuttle receiver. This occurred only when the chase planes were almost 
in-front (position) of the receiver. 

The rain model discussed in Section 7 was installed at the bottom of the 
water tank, and compressed air at predetermined settings was passed through 
the air stones of the rain generating set in order to create a volume of ran- 
dom size bubbles simulating rain. Theoretical and practical justification 
for such a model is given in the next section. 

The Scientific Atlanta positioning system is synchronized with the Sci- 
entific Atlanta Receive and Record System located at one end of the water 
tank by transmit and receive synchros. Thus the recording chart paper can be 
either manually or automatically controlled from the main controller. The 
scan speeds can be set by the setting of the chart speed controller, and the 
calibration of these speeds in degrees per second is given in Fig. 6.4. It 
was found necessary to operate at approximately 1.75 degrees per second in 
order to capture all the details and yet not be affected by the mechanical 
lag of the recording system. 

The experimental data was taken by setting up the system for a given 
range, and the corresponding receiver tilt angle and elevation, under no rain 
conditions, and then repeating it under two rain settings R^ and R^ as shown 
in the rain runs for various ranges in the Appendix. A second receive trans- 
ducer was located to the side of the one simulating the Orbiter receiver such 
a way as to be completely outside its main beam. This unit is called the 
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reference unit, and its position is kept fixed during the entire set of ex- 
perimental runs. The reference transducer is electronically switched in after 
recording the major lobe as well as the following first null of the Orbiter 
receiver beam. The basic purpose of this receiver is to locate the position 
of -4 db beam center of the Orbiter receiver beam from a known and fixed 
reference point, thus enabling us to calculate the shift in the shuttle beam 
due to any multipath distortion may be called a location error for both azi- 
muthal and elevation scan results. A typical calculation of the same is shown 
in Fig. 6.5. 

A summary of all data analysis based on the above criterion is given in 
Table 6.1. An examination of the above data shows that ultrasonally simulated 
multipath errors in azimuthal and elevation shuttle scan beam locations are 
limited to a maximum of 0.08 degrees at 27,250 feet range for azimuthal scan 
and 0.125 degrees at 32,000 feet range for elevation scan. 

A summary of errors analysis as obtained from simulation is given in 
Table 6,2, It is necessary to discuss some of the nomenclature used xn re- 
porting these results. The "NO-RUNWAY" condition implies that the receiver 
and transmitter are located at exactly the same relevant positions in free 
space as in the case of MSBLS in the presence of runway with its surrounding 
ground area. This data point is often used to determine the ground multipath 
as compared to the free space case. Furthermore the data is recorded in 
Tables 6.1, 6,2 and 6.3 in two basic formats; 

Table 6.1 Actual -4 db beam points for both reference and MSBLS simu- 
late receivers. 

Table 6.2 Shift of -4 db beam center of the simulated MSBLS receiver 
with reference to that of the reference beam center (AR) . 

Table 6.3 Actual variation of the -4 db beam width of the simulated 
MSBLS receiver. 
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In view of the above, one has the option to examine the data from different 
viewpoint of the user. 
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TABLE 6.1 SWmRY OF SIMULATED MSBLS DATA AND ANALYSIS 


RUN 

s 

TYPE 

RECEIVER 

' REFERENCE 



DA 

a-/, 

“+4 


; a^.4 

80209 

DIST 

AZ 

SCANS 


iV - 

^ NO R 

AIN 


001 

26.3" 

l^9i 


IBIHH 

mm 

18.5 

I 60.5 

1 A2 

002 

35.3 

KRii 



BM 

1 

15 

1 61 


003 

51.8 

8,5 

18.5 

160.5 

169 

13.5 

16A.75 

IBHBB 

0O»( 

68.8 

0 

1 2 

1 50.5 

' 161 

6 

'155,75. . 


005 

EBBMi 

} k .5 

28 

1 69 

177 

1 21.25 

: 173 

1 151.75 

006 

102.2 

-3 

A 

i IA6.5 

^ 155 

! “2.5 

' 350.75 

I 153.25 

007 

120 


0 

1 138.5 

1 A9 

-8 

^ IA3.75 

1 151.75 

8011 7 


M 

SCANS 


' -R’.IN 



j 

021 

mm 

■»: 

16 

1A7 

i 160 

10.25 

' 153.5 

IA3.25 

022 


-7 

k .5 

1 A2 

153 

-1.25 

’ 1A7.5 

RBBWB 

023 


BHI 

19 

158.5 

1 

1 171 

13.25 

16A.75 

IBB 

02 k 

jggai 

-3.5 

BSRi 

lAg 


mmm 

155.5 

153.5 


Mil 

-3 

HBi 

1A9 



.155.5 

153.25 


.82 

-3 

7.5 

^k3 

! 162 

2.25 

155.5 

- 153.25 

025 


-1 

1 0 

152 

1 62 

A, 5 

1 157 

i 152.5 


Rl 

-1.5 

1 0 

IBnH 

161,5 

A. 25 

' 156.5 

152.25 


87 

L- i ■ 5 

' 10.5 

152 

162 

A. 5 

1 157 

152.5 

026 


wm 


170 

wmm 

20.75 

! I7A.5 

153.75 



13.5 


168,5 


20 

i 173.75 

1 153.75 


R2 

12.5 

mHi 

MEmam 


19.5 

173.75 

1 15A.2 5 

027 

' ^9.^nr 

19 

3 A 

173 

182 

26.5 

' 177.5 

! 1 51 


R, 

18 

35 

17.3 . 

183 

26.5 

msEmmmmm 

I 151.5 


mamm 

- 

- 

173.5 

TP 

- 

178.75 

- 

SoTTTo 



EL 5 


/AY - 

RAIN 



028 

mm 

5.5 

11.5 

32,5 

38 

8.5 

35.25 

26.75 

029 

mm 

0 

5.5 

55 

61 1 

2.75 

58 

55.25 

030 


-It 

1.5 

60 

65 

-1.25 

62.5 

63.75 

; j 

Ri 

-A .5 

1.5 

59 

65 

-1.5 

62 

63.5 

1 

«7 

' 5 

1.5 

59 

65 

-1.5 

62 

63.5 

031 

mam 

m 

BRH 

57,5 

63 

-7. 75 _ . J 

60.25 

68 


ilMHi 

‘[mil 

mm 

57.5 

63 

-7.75 

60.25 

68 


P2 



57.5 

63 

-7.5 

60.25 

67.75 

032 

iicgaa 

-3.5 

3.5 

68 

7 k 

0 

71 

71 




3 

67.5 

73 ! 

-0.5 

70.25 

70.75 



-k 

3 

67.5 

73 

-0.25 

70.25 

70.5 

553 


"875“ 

li 



«7.5 

11.75 

■8 A": 7 5 

73 


Rl 

8.6 

_L5 

_U 

--87 

-1 U.7.5 

_a_A. .. . . .. 

72.. 2 5 
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TABLE 6.2 SUMMARY OF SIMULATED MSBLS MULTIPATH 


CAUSED -4 db BEAM (IN REFERENCE TO THE CORRESPONDING REFERENCE BEAM CENTER) 


Actual 

Range 

Azimuth Scan 



Elevation Scan 


(kilo ft) 

NR 


^2 

NR 


^2 

7.0 

+0.209 



Ref 



9.4 

+0.455 



Ref 



13.8 

+0.043 



Ref 



18.35 

+0.625 

-0,042 

0.000 

Ref 

-0.042 

-0.042 

22.85 

+0.125 

-0.042 

0.000 

Ref 

-0,000 

-0.042 

27.25 

+0.083 

0.00 

0.083 

Ref 

-0.042 

-0.083 

32,0 

-0.125 

+0.081 


Ref 

-0.125 

-0.125 


NR = NO RAIN CONDITION 
R^ = RAIN RATE ONE 
R^ = RAIN RATE TWO 
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TABLE 6.3 SUMMARY OF SIMULATED MSBLS RECEIVER -4 db 


BEAM WIDTH (x6) 


Actual Range 
(kilo ft) 

NR 

Azimuth Scan 
^1 

^2 

NR 

Elevation Scan 

^2 

7.0 

12.5 

A 

A 

AA 

A 

A 

9.4 

11.5 

A 

A 

6 

A 

A 

13.8 

11.5 

A 

A 

5.5 

A 

A 

18,35 

11.0 

10.5 

10.5 

5.5 

6.0 

6.0 

22.85 

11.0 

11.5 

12.0 

6.5 

6.5 

7.0 

27,25 

13.5 

13.0 

14.0 

7.0 

7.0 

7.5 

32.0 

15.0 

17.0 

aaa 

6.5 

6.5 

7.0 


* - Rain could not be simulated because of air bubble accumulation 
on model runway, 

** - Physically not possible to run because of positional problem 
of equipment. 

AAA - Operational problems - data is questionable. 
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7. ULTRASONIC SIMULATION OF RAIN 


Multipath degradation of azimuth and elevation scan signals for ground- 
based transmitters and airborne receivers such as MSBLS are expected to get 
worse for light and heavy rainfall conditions as opposed to no-rain conditions. 
The maxn reasons for this effect are the different phase fluctuations of the 
direct, and indirect path signals because of 

a) the intervening rain volume Instead of free space, 

b) the nonuniform rain distribution along both paths, 

c) the propagation through rain volume at different angles (measured 
from the vertical) for both paths, 

d) the additional random fluctuations of the amplitude and phase angle 
of the forward reflection coefficient from the ground points due 

to the time-varying ground surface disturbed by raindrops during 
rainfall, 

e) the variation of the average attenuation and composite phase shift 
for a uniform rain volume path length. 

It is often considered sufficient, in case of line-of-sight propagation paths, 
to assume the path attenuation, for uniform rain for a given propagation fre- 
quency to be: 

total db attenuation = No rain attenuation + rain caused attenuation 

= 20 log (X/4irr) + KR (7.1) 

where K = constant db/n. miles, and it varies with rainfall rates and 
the operating frequency 
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R = range - n. miles 


For instance, at a frequency of 15 gigahertz, one obtains the value of K as 
1.67 db /nautical mile for 0.63 inches per hour of rainfall. This is quite 
satisfactory for line— of— sight communication paths where one is interested 
in long-time averaging of signals, but MSBLS requires accuracies of 0.05 
degree, (equal to one standard deviation) and MSBLS updates its average 
positional readout averaging period approximately every 200 milli-seconds. 

It is quite obvious that long-time average attenuation or phase changes would 
not be very helpful in determining iKBLS multipath degradation in the pre- 
sence of rain, let alone the multipath amplitude and phase effects of the 
rain effected ground as is the case for MSBLS. 

It is also very pertinent to note that in all theoretical work in the 
area of rainfall attenuation and depolarization (Warner et al, 1976, Mink, 
1976) , rain drop shapes and their orientations models have been assumed for 
uniform rain drop size distribution and rain rate along, the path. This 
assumption of uniform rain rate, and the use of average rain rate is ques- 
tionable in the case of Kennedy Space Center environment which has experienced 
a very short duration squall of 11,52 inches /hour rain for 0.3 minutes< 

(Gullick 1976). It is believed that during such extreme rainfall periods, 
one may encounter the worst case MSBLS multipath degradation. In order to 
predict the degradation, the rainfall is simulated by air bubble volume in 
water at room temperature, in which air bubble size and shapes are randomly 
distributed just as is the case with the rain drop sizes in the natural en- 
vironment, First the theory of rain modeling is discussed and then are 
given experimental results of rainfall attenuation for different rain rates 
in order to have a base for comparison with full scale results. Furthermore, 
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the beam shape distortion caused by ground multipath for MSBLS, although 
smoothed out by faster scan rates, leaves an unsymmetrical beam which is 
considerably different for that obtained for the case of the line-o£-sight 
propagation in rain volume. 

RAIN DROP MODELING: AIR BUBBLE SIZE DETEEMINATION 

The modeling of rain drops is based on simulation of their radar cross- 
section and transmission coefficient of rain volume. The rain drops are 
simulated by air bubbles in a water tank, and acoustic sources with appro- 
priate bearawidths are used to simulate electromagnetic sources of MSBLS. 

The air bubbles are formed by the flow of pressurized air through a set of 
fish stones located at the bottom of the tank. Let us first examine the 
details of the radar reflectivity of a static rain model. 

The velocity of the falling speed of water for a rain drop of radius a 
cm is (Gunn and Kinzner, 1944) 

V(a) = 950 a-exp(-a/0.875)^*^} cm/sec (7.1) 

Eq, (7.1) may be approximated by the following Eq. for the range 0.5 < a < 
0.3 cm. 

V(a) = 950 {l~exp(-v a)} cm/sec (7.2) 

where y = 11.2 cm ^ 

The radar extinction cross-section of a single drop (Ryde and Ryde, 1944) 
can be approximated by 

= 120a^'^^ cm^ (7.3) 
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and the extinction cross-section for the range .03 < a < .25 cm is given in 
Table 7.1 for X = 3.2 cm. 


The attenuation Y due to the rain drops is given below. 


Y = 1 N(a) o . (a) da (7.4) 

Jo 

where N(a) is the distribution function of the size of the rain drop and 
a (a) is defined in Eq. (7.3). N(a) is assumed to follow the distribution 
function given by Khrgian and Mazin (1952) or 


N(a) = CaV^^ 


C = 54 cm 


(7.5) 


where g is a parameter (cm ^) significant of the rain intensity, R, which is 
given as follows: 


R = 


r N(a, 


Jo 


4ira^ 

3 


V(a) da 


(7.6) 


Substituting V(a) and N(a) from Eqs. (7.2) and 7,5) in Eq. (7.6) one obtains 


R(?) = 7.75 X 10^ |~ - (7.7) 

The values of R and y for various values of C are given in Table 7.2 (Kwan, 
1968) for X = 3.2 cm. 


The backscattered received power Pr is given by the following: 
^-mycT (L^-L^) g2^2 


p ^ p ^ 

Y t Li 


(4tt)‘ 


f (0,(#i)dij) d0 

Jo 


aN(a) da 


(7.8) 


where G = antenna gain 
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T = the pulse length 
m = an integer » 1 
(L2“L^) = thickness of the rain layer 
X = wavelength 

= transmitted power 

O 

C = velocity of light 3x10 m/sec 
$,0 = antenna beamwidths 
- antenna beam 


Substituting (Battaan, 1959) the value of G given below 


G = 


$9 

//f^(0,(|))d^de 


1.76 


where 


2 2 

- exp - (-^ + -^) 4 Ln 2 
0 $ 


(ca “> W 81 

N(a) a da = ^ 1 

n=-co (g-inw)" 


, e -1 2 
4 I c 


t|) = 4 -ft (27T/X) l-^^l 

c 


W = the Fourier coefficients (Stephens, 1960) in 

■ I """ 

n=-“ 

= complex permitivity of rain volume 


in Eq. (7.8), one obtains: 


P 

r 



1.76 

0$ 


“mycT 


h2-L 


i tc I 


n=-co 


W 81 
n 

(5-inw)^ 


(7,9) 

(7.10) 

(7.11) 


(7.12) 
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The Mie Backscatter cross-section is given by: 


CT 


mie 


0 . 

mie 

a 

ray 


0 

ray 



I 


W 

n 


^jnw 


(7.13) 


0 . 
lUXS 

where — — is tabulated by Stephen (1960) , for a propagation wavelength of 
ray 

3.2 cms. in Table 7.3, and one can write the following results: 


0 

ray 


0.95 



= 1.83 X 10^ 


6 

a 


for X - 3.2 cm 


(7.14) 


The extinction cross-section of axr bubbles for an acoustic source in 
water is (Urick, 1968) given below: 


4it 5/K a 

r 

^e ~ 2 2 

[(f^/f) - 13+6 

where f^ = resonant frequency -Hz 

6 = damping constant 

a = radius of bubble -cm 

K = 27Tf /C„ 
r Y 0 

Cq ^ velocity of sound in water 


(7.15) 


The resonant frequency for a bubble is given as 


f 

r 


326 

a 


[1 + .03d] 


1/2 


(7.16) 


where d the depth of the water in feet, and a is measured in centimeters, 

Tor this rain simulation, d is assumed to be 3 feet and hence f may be 

Y 

approximated by 326/a and K^a by .0186. In order to simulate the study of 
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radar cross-section due to rain drops, the size of air bubbles must be such 
that it gives the same scattering cross-section for the acoustic source. 

The radius of the air bubble a is given by the following equations; 


6 4 2 

a + + Mj^a + = 0 

where 

= -0^(1 + 6^)/4ir 
M = 2 X C326)^ a^Airf^ 

1 D 

= -cr (326)^/4irf^ 

U e 


(7.17) 


(7.18) 

(7.19) 

(7.20) 


The coefficients and are a less than unity at the frequencies of oper- 
ation, viz, .01 mhz < f < 15 mhz, and therefore (7.16) becomes: 

a^ + M^a^ = 0 (7.21) 

which has only real roots. An analytical solution of this Equation yields 
the bubbles radius given below: 


a = /a, /4ir for a X (7.22) 

b 

The size of the air bubbles for the frequency of 2.25 mhz are given in Table 
7.4. 

Furthermore one may use the coefficient of attenuation K in db/cm de- 

e 

fined by the equation given below. 


K = 4.34 n a (7.23) 

e e 

where 

~ 3 

n = average number of bubbles per cm 
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2 

ff = extinction cross-section of bubbles in cm 
6 

It is most pertinent to note, that one employs the probability distribution of 

diameter of bubbles in the case of underwater bubbles in the same manner as that 

for raindrops in the case of radar. The expression (7.23) may be modified to 

include the total volume of an V(a) contributed by bubbles in the beam in water, 
to obtain 

K = 1.504 X 10^ v(a) (7.24) 

This expression is analogous to total rain drop contributed volume of water 
in air for the case of radar. 

CALIBRATION OF RAIN SIMULATOR 

Since radio and radar frequency attenuation by rain volume has been 
usually reported in the literature in the form of average db/km attenuation. 

For instance Livingston (1970) reports an attenuation of 1.243 db/km for 
15 mm/hour raxn rate at 16 GHz, whereas Kerr (1951) summarizes rain atten- 
uation measurements at the same frequency and for the same rain rate to be 
1.65 db/cm. Furthermore McCormick et al. (1976) report that attenuation 
rates of 1 to 15 db/km at 16.5 GHz measured over a long period of observation 
varying from 1972-1975, and further show the importance of canting angle. 

The basic reasons for such a broad range of variations are: 

1) the rain volume attenuation measurements are averaged over many 
minutes and in fact up to an hour or so at many times 

2) the rain volume is assumed to have a uniform ram rate 

3) the extent of rain free axr volume in the radar path is not known 
to the experimenters 
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4) there is no ground multipath, as in the case of MSBLS, which 
could further broaden the range of attenuation. 

In view of the above, it was decided to use an average attenuation figure 
from Livingston (1970) for 15GHz frequency, to calibrate the rain simulator. 
Airstones spanning a simulated range of 28 inches (corresponding to a full 
scale MSBLS range of 2.28 kilometers) were used to generate continuous rain 
volume and the transmitter/receiver were set up in their usual location in 
the MSBLS geometry so that the rain caused fading statistics would be more 
meaningful than that affecting a radar beam traversing a rain volume at 
right angle to the path of the rain drops. The theoretically verified ex- 
perimental rain attenuation data at 15 GHz used in this calibration is 
given below: 


Rain Rate mm/hr (IN /HR) Attenuation db/Km 


4.0 

(0.16) 

0.30 

10.0 

(0.39) 

1.15 

20 

(0.79) 

2.4 

30 

(1.18) 

3.6 

40 

(1.57) 

4.5 

50 

(1.97) 

5.3 


Based on the above reference data, the simulated rain rates and R 2 were 
found to be equivalent to 0.31 "(7.8 inm/hr) and 0.61" (11.5 mm/hr) inches 
per hour respectively, since attenuation caused by rain only of 2 and 4 db 
were found to occur for a simulated 28" (equivalent to 2.276 Kn full scale) 
long rain volume along the path of the simulated radar beam. The canting 
angle of the simulated radar beam was taken into account in calculating the 
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equivalent rain rates. 


FADING DEPTH 

The multipath caused fading depths are measured in db excursions from 
the average signal level. Both experimental and theoretical rain caused 
fading data for 15 GHz signals (Goldhirsh 1975) show that, for a thunderstorm 
or cloudburst (approximately 100 mm/hr or 3. 94" /hr rain rate) , the probability 
P(d>D) of the fade depths exceeding a given fade depth D for MSBLS radar 
beam tilts may be as given below; along with projected values for 1"/HR rain: 

100 mm/HR 25.4 mm/HR 


20.0 

2.56 

0.56 

15.0 

5.12 

2,58 

6.0 

10.24 

5.18 

3.5 

15.39 

10.3 

2.5 

20.48 

15.4 

1.8 

25.6 

20.5 

1.5 

30,78 

25.7 


The fade depths for simulated rain rates of R^ and R 2 were found to be 3.13 
and 7.88 db respectively, which were exceeded, approximately 6 and 3 percent 
of the time. These fading statistics fall within the expected ranges calcu- 
lated from data and discussion by Goldhirsh (1975) and Crane (1975). 

MSBLS RAIN MULTIPATH SIMULATED RESULTS 

Two simulated uniform rain rates and R 2 , described above, were used 
to collect data on MSBLS receiver beam shape in order to study the multipath 
effect on its -4 db width. Only 28" (2.276 Km in full scale case) of the 
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total horizontal path length contained rain, and this extent of simulated 
rain volume was approximately 23% of the maximum simulated horizontal range 
( 120 "). 

In the case of MSBLS ultrasonic simulation, the radar signal back-and 
forward-scattered by rain volume is appropriately reproduced in simulated 
rain in the form of air bubbles in water. For instance, MSBLS guidance sig- 
nal at approximately 15.55 GHz (1.93 cm, wavelength) is simulated by 2.25 mhz 
(0,667 mm wavelength) ultrasonic azimuth and elevation scan signals with 
appropriate waveform and beamwidths. The time varying wind as well as rain 
intensity achieved by controlling the compressed air flow used for generating 
bubbles and simulating wind gusts by a jet of water blown into this volume 
of air bubbles were not Incorporated because of severe fading obtained with- 
out these. 

The minus four db receiver bearawidths for the simulated MSBLS case are 
listed in Table 7.5. The effect of the presence. of rain apparently is to 
widen the said beamwidth in general by up to 2 degrees, whereas in few cases 
a reduction of up to 0.5 degree takes place. This is significant since these 
changes correspond to up to +13.3% increase&of -4.5% decrease respectively 
of the no rain condition case values. In the case of one inch rate, one 
would expect these changes in minus 4 db bearawidths to be more pronounced 
and of the order of up to +23% and probably no reductions. The ultrasonic 
simulation provided very extensive fading data discussed above. 
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TABLE 7.1 BUBBLE SIZE AND ITS EXTINCTION CROSS SECTION 


ext^ 2 

Bubble / Drop Size (cm) Log j^^ e cm 


.03 

6.5 

X 

io “^ 

.04 

2.5 

X 

10 "^ 

.05 

7 

X 

10 -^ 

.06 

1.5 

X 

lO -"^ 

.08 

5 

X 

10 "^ 

.10 

1.5 

X 

10-3 

.15 

9 

X 

10-3 

.20 

3 

X 

10-3 

.25 

8 

X 

IQ -3 


TABLE 7 : 2 RAIN INTENSITY (R) AND ATTENUATION (db/m) 


R 

mm/hr- ^ 


(m~^) 


10 Log _e ^ 
db /meter 


170 

40 

1.1 

X 

10 -=> 

4,75 

X 

10-3 

100 

43.7 

.56 

X 

10-3 

2.43 

X 

10-3 

43 

50 

. 2025 X 

10-3 

0.88 

X 

IQ- 3 

13 

60 

.51 

X 

1 

o 

2.21 

X 

IQ-'^ 

4.66 

70 

,16 

X 

IQ-'^ 

.694 

X 

10 -^ 

1.93 

80 

.58 

X 

10"3 

2.51 

X 

10“3 

0.88 

90 

.236 

X 

10~3 

1.02 

X 

IQ-3 
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TABLE 7.3 RAIN-DROP BACKSCATTERING CROSS SECTION 


a 

ray 


(Backscattered Cross Section) 


.05 

1 

2.75 

X 

10-6 

.275 

X 

10-3 

.075 

0.9 

.32 

X 

10"'^ 

.28 

X 

IQ-''^ 

.1 

0.85 

.183 

X 

10-3 

.16 

X 

10“ 3 

.125 

0.75 

.7 

X 

10-3 

.53 

X 

10-3 

.15 

0.80 

.21 

X 

io“3 

.17 

X 

10"3 

.175 

1.35 

.53 

X 

10“^ 

.72 

X 

IQ- 3 

.2 

2.35 

.116 

X 

10"^ 

.27 

X 

IQ-^ 

.225 

2.4 

.24 

X 

IQ-^ 

.0576 



.25 

2.4 

.44 

X 

IQ-^ 

.105 



.30 

1.8 

.96 

X 

10~3 

.173 



.35 

1.65 

.335 



.552 



.40 

1.4 

.74 



1.04 



.45 

1.1 

1.484 



1.63 



.50 

1 

2.97 


- 

2.97 




TABLE 7.4 AIR BUBBLE RADIUS FOR DIFFERENT BACKSCATTERED 
CROSS-SECTIONS AT 2.25 MHZ 


Desired Backscattered Crosssection Bubble Size - xnm 


12.57 1 

50.28 2 

113.13 3 

201.12 4 

314.25 . 5 

452.52 6 

615.93 7 

804.48 8 

1018.17 9 

1257 10 

2828 15 

5028 20 
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TABLE 7.5 MINUS FOUR DB RECEIVER BEAMWIDTH (DEGREES) OF 


ULTRAS ONICALLY SIMULATED MSBLS 


Simulated Range (in) 
Conditions 

26.3 

35.3 

51.8 

68,8 

85.7 

102.2 

120 

AZ - No Rain 

11.5 

11.5 

11.5 

11.0 

11.0 

13.5 

15.0 

AZ - Rain R^ 

— 

— 

— 

10.5 

11.5 

13.0 

17.0 

AZ - Rain R^ 

— 

— 

— 

10.5 

12.0 

14.0 

— 

AZ - Rain R^ 

12.0 

12,2 

12.5 

13.0 

13.5 

15.0 

18.5 

EL - No Rain 

6.0 

5.5 

5.5 

5.5 

6.5 

7.0 

6.5 

EL - Rain 

— 

— 

— 

6.0 

6.5 

7.0 

6.5 

EL - Rain R^ 

— 

— 

— 

6.0 

7.0 

7.5 

7.0 

EL - Rain R^ 

6.1 

6.0 

6.3 

6.5 

7.5 

8.0 

7.5 


Notation 

= 0.31"/HR (7.8 mm/HR), R2 = 0.61" /HR (11.5 mm/HR) 

R^ = 1"/HR (25.4 mm/HR) Projected values of beainwidth changes. 

— « No rain data available because of accumulation of air bubbles on 
simulated runway surface. 
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FIG. 7.1 RAIN GENERATION AIR STONES SETUP 
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8. Conclusions 


MSBLS antenna coverage, landing trajectory and the ground plan contain- 
ing the runway and its surrounding area were ultrasonically simulated. Fixed 
points were selected on the- trajectory for collecting MSBLS simulated fading 
data for no rain, and two rain rate conditions of 0.31" and 0.61" per hour. 
Rain was simulated by air bubbles of random radii generated by air stones in 
water . 

Rain caused minus 4 db beamwidth fading for 0.31"/HR simulated rain (R^) 
was found to vary from -4.5 to 13.3%, and 0 to 9.1% for Az and El - beams 
respectively. The corresponding range from -4,5 to 16.7%, and 7,1 to 9.1% 
respectively were obtained for 0.61" /HR rain (^- 2 ^ • Similar projected values 
for 1.0" /HR (R^) were found to be 4.4 to 23.3% and 1,7 to 14,2% for AZ and 
EL beams respectively. 

The fade depths with a probability of occurrence of 15% for R^, R^ and 
R^ rain rates were obtained as 0.5, 0.9 and 2.56 db respectively, whereas 
those for 6% occurrence these values become 3.6, 4.5, and 5.18 respectively. 
It is noteworthy to add that fade depths of 24.1, 25, and 25,7 db may occur 
1.5% of the time for R^, R^ and R^ respectively. 

Cloudburst, common to Cape Kennedy, defined as 100 mm/HR (3. 94" /HR) 
may cause a fade depth of 2,56 db twenty percent of the time and 5.12 db 
fifteen percent of the time. The corresponding minus 4 db receiver beam- 
widths are expected to fluctuate far more severely for this case, such as 
30-50% of the rain beamwidth values. 

The fade depths as well as fluctuation of the simulated minus 4 db 
receiver beamwidth don’t seem to follow a systematic range variation because 
of the nonlinear relationship of the illuminated ground area. Another sig- 
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nificant result of this simulation is that the simulated receiver beam shape 
is considerably more distorted by rain than indicated by the minus 4 db beam- 
width variation. This will effect the MSBLS receiver output subject to the 
analysis process, and its sensitivity, since a jagged beam will display 
sudden changes in slope from one sample pulse to the next, and the orbiter 
system may just keep dropping off such sample values and if such severe 
environmental conditions continued to recur in a complete cycle, the whole 
cycle may result in a complete blanking of the guidance data for the said 
cycle. Of course, if the MSBLS receiver smooths out the input data, then 
it may overcome such major fluctuations in the signal amplitude caused 

by degradation of the recexver beam shape. 

The actxjal percent path length containing no rain can make a considerable 
difference in the validity of the standard rain attenuation formula and 
associated average rain attenuation rates. In real time fast systems such 
as MSBLS, there is considerable fluctuation of the signal, and it is anti- 
cipated that MSBLS sampled data points taken on a beam would experience all 
the random fluctuation of the rain attenuation as well as those caused by 
the phase and amplitude fluctuations caused by ground multipath as opposed 
to that experienced by a direct path. 

The effect of chase aircraft as a multipath degradation defect was 
found to be negligible, in fact the recorder could not record any noticeable 
difference in the beamshape as compared to the case without chase aircraft 
due to the narrowness of the azimuth elevation fan beams. 

Finally the results obtained in the ultrasonic simulation at the Uni- 
versity of Houston - Wave Propagat^ion Lab are believed to a reasonably pre- 
diction of the full scale results. 

It is recommended that a detailed crosscorrelation field data, and UH 
simulated data should be performed. 
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10. APPENDIX 



LIST 

OF RUNS AND RUN 

DESCRIPTION 

RUN NO. AZ/el 


rain (R] , R2) 

FREESPACE (FS) 

RANGE (k-ft.) 

NO RAIN (NR) 

RUNWAY (RW) 

80209001 -AZ 

7.0 

NR 

FS 

9002-AZ 

9.^ 

NR 

FS 

9003-AZ 

13.8 

NR 

FS 

9004-AZ 

18.35 

NR 

FS 

9005-AZ 

22.85 

NR 

FS 

9006-AZ 

27.25 

NR 

FS 

9007-AZ 

32.00 

NR 

FS 

801 17021-AZ 

7.00 

NR 

RW 

1022-M 

9. A 

NR 

RW 

7023-AZ 

13. e 

NR 

RW 

7024-AZ 

18.35 

NR,R, ,R2 
NR,R^ ,R2 

RW 

7025-AZ 

22.85 

RW 

7026-AZ 

27.25 

NR,Rj ,R2 

RW 

7027-AZ 

32.00 

NR,R| R^ 

RW 

7028-EL 

9.4 

NR 

RW 

7029-EL 

13.8 

NR 

RW 

7030-EL 

18.35 

NR, R^,R2 
NR,R, ,R2 
NR,R^,R2 
NR,RpR2 

RW 

7031-EL 

22.85 

RW 

7032-EL 

27.25 

RW 

7033-EL 

32.00 

RW 
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